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ABSTRACT 

The  Naval  Research  Laboratory  (NRL)  has  developed  a  global,  relocatable,  tide/surge  forecast  system 
called  PCTides.  This  system  was  designed  in  response  to  a  U.S.  Navy  requirement  to  rapidly  produce  tidal 
predictions  anywhere  in  the  world.  The  system  is  composed  of  a  two-dimensional  barotropic  ocean  model 
driven  by  tidal  forcing  only  or  in  conjunction  with  surface  wind  and  pressure  forcing.  PCTides  is  unique  in 
its  ability  to  forecast  tidal  parameters  for  a  user-specified  latitude/longitude  domain  easily  and  quickly,  and 
is  especially  useful  in  areas  where  observations  are  nonexistent.  PCTides  provides  short-term  (daily  to 
weekly)  predictions  of  water-level  elevation  and  depth-averaged  ocean  currents.  The  system  has  been  tested 
in  numerous  regions  and  validated  against  observations  collected  in  conjunction  with  several  navy  exercises. 


1.  Introduction 

The  ebb  and  flow  of  tidal  waters  are  of  importance 
for  a  wide  range  of  military  applications  including  mis¬ 
sion  planning,  amphibious  landings  (e.g.,  the  D-Day  in¬ 
vasions  at  Normandy,  France;  and  Inchon,  South  Ko¬ 
rea),  ship  routing,  navigation  through  harbors,  bathym¬ 
etry  processing  (i.e.,  the  removal  of  tidal  effects  from 
sounding  data),  and  the  operation  of  unmanned  under¬ 
water  vehicles  in  shallow  water.  Tidal  currents  can  af¬ 
fect  the  positioning  and/or  movement  of  subsurface 
mines  and  the  performance  of  U.S.  Navy  sea,  air,  or 
land  (SEAL)  swimmers. 

The  U.S.  Navy’s  globally  relocatable  tidal  prediction 
model,  PCTides,  was  developed  to  fill  a  void  in  the 
navy’s  global  tide  forecasting  capability.  Previously, 
tidal  forecasts  available  to  the  navy  were  confined  to 
coastal  locations  where  water-level  data  were  available. 
This  restricted  tidal  forecasts  to  small  areas  between 
tidal  stations.  In  addition,  these  forecasts  did  not  in¬ 
clude  the  effects  of  wind  (surge),  which  can  play  a  sub¬ 
stantial  role  in  water-level  prediction. 

PCTides  was  developed  as  a  global  tide-forecasting 
system  that  can  be  set  up  rapidly  for  any  user-specified 
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location(s)  to  predict  the  water-level  elevation  and 
depth-averaged  ocean  currents.  In  addition,  PCTides 
can  assimilate  the  local  meteorological  effects  of  wind 
and  surface  pressure  gradients,  which  typically  play  a 
role  in  water-level  prediction.  PCTides  was  specifically 
designed  to  produce  daily  to  weekly  forecasts  of  tidal 
conditions;  however,  longer  predictions  (e.g.,  months) 
can  be  made  for  military  mission  planning  as  well  as  for 
estimating  optimal  times  for  deployment  of  instrumen¬ 
tation,  etc.  The  global,  rapidly  relocatable  tide/surge 
forecast  system  is  implemented  on  Windows  operating 
systems  (e.g.,  XP)  running  in  a  DOS  mode  as  well  as  on 
UNIX  and  LINUX  platforms. 

Section  2  of  this  paper  describes  the  PCTides  system, 
and  section  3  presents  a  discussion  of  validation  tests 
performed  at  various  worldwide  locations.  A  summary 
and  conclusions  are  presented  in  section  4. 


2.  The  PCTides  system 

The  PCTides  system  has  a  two-dimensional  (2D) 
barotropic  ocean  model  as  its  core  component  (Fig.  1 ). 
PCTides  uses  this  barotropic  ocean  model  to  predict 
both  depth-averaged  currents  and  sea  level  heights  on 
or  near  continental  shelves  (Preller  et  al.  2002).  It  con¬ 
tains  a  wetting  and  drying  algorithm  for  the  simulation 
of  coastal  flooding  due  to  tides  or  storm  surge.  The 
model  is  driven  by  tidal  forcing  and,  if  desired,  atmo¬ 
spheric  surface  wind  and  pressure. 
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Fig.  1.  The  PCTides  system  components. 


a.  The  PCTides  Ocean  Model 

PCTides  uses  the  following  shallow-water  equations 
to  define  the  depth-averaged  current  [Eqs.  (l)-(2)]  and 
the  sea  surface  elevation  [Eq.  (3)]: 
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chosen  map  projection  of  the  model  grid);  #  is  the  ac¬ 
celeration  due  to  gravity;  P  is  the  atmospheric  surface 
pressure;  pM  is  the  water  density;  v  is  the  coefficient  of 
lateral  eddy  diffusion;  Thx  and  rhy  are  the  bottom  fric¬ 
tional  stress;  tsx  and  rvv  are  the  surface  wind  stress  in  the 
x  and  y  directions;  and  Sxxf  Svv,  Sxy,  and  Svx  represent 
surface  wave  radiation  stresses. 

The  surface  wind  stress  components  are  computed 
using  the  quadratic  relationship: 

Tsx  =  CnP„  I  “  ,  I  U,r  T.v.v  =  CnP„  I  U  ,  I  v,r  <4» 


where  U  and  V  are  the  depth-averaged  currents  in  the 
x  and  y  directions,  respectively;  H  is  the  total  depth;  £  is 
the  sea  surface  elevation; /is  the  Coriolis  parameter;  m 
is  the  map  factor  (a  scaling  factor  dependent  on  the 


where  luj  =  ( u l  +  i^)(1/2>,  ua  and  va  are  the  horizontal 
components  of  wind  velocity  at  anemometer  height,  p(l 
is  the  density  of  air,  and  CD  is  the  drag  coefficient  based 
on  Smith  and  Banke  (1975),  which  is  expressed  as  fol¬ 
lows: 

CD  =  [0.63  +  0.066| uj]  x  10-3  |uj  <  25  m  s1 
CD  =  [2.28  +  0.033( | u J  -  25)]  x  10"3|uJ  >  25  m  s”1. 

(5) 

The  bottom  stress  is  represented  by  a  Manning’s  n 
depth-dependent  friction  relation  (Signell  and  Butman 
1992): 
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Fig.  2.  Bathymetry  (m)  of  the  Yellow  Sea.  Squares  indicate  the  IHO  station  locations  and  circles 
indicate  the  location  of  the  ADCPs  during  September  1995. 


8^ 

Thx  P"  (H  + 


Tby  =  Pw 


(H+Q 


1/3 


(U2  +  V2)U2U, 


(U2  +  V2)v2V,  (6) 


where  n  has  the  value  0.03.  This  formulation  ensures 
that  the  drag  coefficient  increases  with  decreasing  wa¬ 
ter  depth  and  is  applied  to  water  depths  greater  than  1 
m.  In  extremely  shallow  water  and  over  land  points  that 
become  inundated,  the  bottom  drag  coefficients  can  be 
specified  at  each  grid  point  according  to  the  terrain 
type. 


Equations  (l)-(3)  are  solved  on  an  Arakawa  C  grid 
(Messinger  and  Arakawa  1976)  using  a  three-level, 
split-explicit  finite-difference  scheme  as  described  in 
Hubbert  et  al.  (1990).  The  continuity  equation  and  the 
gravity  wave  and  Coriolis  terms  in  the  momentum 
equations  are  solved  with  the  shortest  time  step  using 
the  forward-backward  method.  The  nonlinear  advec- 
tive  terms  are  solved  with  an  intermediate  time  step 
using  the  two-time-level  method  of  Miller  and  Pearce 
(1974).  The  surface  wind  stress,  bottom  friction  stress, 
and  atmospheric  pressure  terms  are  solved  with  the 
longest  time  step  using  a  backward-implicit  method. 
This  split-explicit  approach  is  very  efficient  for  oceano- 
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Station  1  (33.2700  N,  125.2510  E) 


Station  2  (33.7500  N,  125.0000  E) 


Station  3  (35.9680  N,  124.0330  E) 


Fig.  3.  Yellow  Sea  PCTides  water-level  elevation  forecasts  (red;  m)  with  data  assimilation  vs  the 
ADCP  observations  (black)  at  stations  1-4  for  the  period  1-30  Sep  1995. 


graphic  models  with  a  free  surface  because  of  the  large 
disparity  between  advective  and  gravity  wave  phase 
speeds  in  deep  water. 

Surface  meteorological  forcing  is  applied  via  the 
wind  stress  and  surface  pressure  gradient  terms  in  Eqs. 
(l)-(3)  at  all  nonland  model  grid  points  in  the  compu¬ 
tational  domain.  Tidal  and  meteorological  forcing  at 
the  lateral  boundaries  is  specified  by  the  incremental 


displacement  of  the  water  surface  due  to  changes  in 
tidal  height  and  atmospheric  barometric  displacement. 
The  lateral  boundary  conditions  are  applied  using  a 
one-way  nesting  technique.  Tidal  elevation  boundary 
values  are  applied  with  decreasing  intensity  from  the 
boundary  to  some  specified  number  of  model  grid 
points  ( nmax  ~  12)  into  the  domain  according  to  the 
following  equation: 
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</>  =  (1  ~  a)4>,,  +  a<f>h.  (7) 

where  <f)p  is  the  model-predicted  value,  <f)b  is  the  pre¬ 
scribed  boundary  value,  and  a  is  varied  according  to  a 
cosine  function  such  that 

a  =  0.5[cos7rtl  -  n/nmax)  +  1],  n  =  1,  nmax.  (8) 

Under  any  prescribed  boundary  forcing,  the  model  is 
initialized  by  setting  velocities  to  zero  and  interpolating 
the  global  Finite  Element  Solutions  version  99  (FES99) 
(Lefevre  et  al.  2002)  elevation  field  to  the  model  grid.  A 
spinup  period  is  used  to  allow  the  effects  of  the  initial 
conditions  and  boundary  forcing  to  propagate  through¬ 
out  the  computational  domain.  Sensitivity  studies  have 
determined  that  the  spinup  time  needed  for  PCTides 
using  wind  and  tidal  forcing  is  24  h  and  is  12  h  when 
wind  forcing  is  not  used. 

b.  PCTides  configuration 

PCTides  was  developed  as  a  stand-alone,  global,  tidal 
prediction  system  and,  as  such,  contains  all  of  the  nec¬ 
essary  global  databases.  These  databases  include 
bathymetry,  tidal  boundary  forcing  data,  and  tidal  sta¬ 
tion  data.  PCTides  uses  a  2-min  global  bathymetry  da¬ 
tabase  developed  by  Naval  Research  Laboratory 
(NRL)  called  NRL  Digital  Bathymetric  database  2 
min  (DBDB2)  (see  http://www7320.nrlssc.navy.mil/ 
DBDB2_WWW/).  This  database  is  derived  from  a 
number  of  sources  including  the  U.S.  Naval  Oceano¬ 
graphic  Office  (NAVOCEANO)  global  dataset 
(DBDBV,  available  online  at  http://gcmd.nasa.gov/ 
records/GCMD_DBDBV.html),  the  Smith  and  Sand- 
well  global  dataset,  the  Data  Assimilation  and  Model 
Evaluation  Experiments  (DAMEE)  North  Atlantic 
data,  the  International  Bathymetric  Chart  of  the  Arctic 
Ocean  (IBCAO)  data,  the  Australian  Bathymetric  and 
Topographic  (ABTG)  data,  and  regional  datasets  from 
the  Gulf  of  Mexico  and  Yellow  Sea.  In  addition,  an 
appropriately  formatted,  user-supplied  bathymetry  can 
be  blended  with  the  default  bathymetry  to  create 
bathymetry  for  highly  resolved  grids. 

The  FES99  global  tidal  database,  which  provides 
tidal  data  for  the  open  boundaries  of  PCTides,  contains 
tidal  harmonics  for  nine  constituents  (M2,  S2,  /V2,  K2a 
Kx,  Ox ,  Px,  Qx,  and  2 /V2).  Although  the  total  tide  is 
defined  by  a  much  larger  set  of  constituents,  these  nine 
provide  a  large  percentage  of  the  total  tide.  The  FES99 
database  has  0.25°  resolution  and  was  derived  by  the 
assimilation  of  approximately  700  tide  gauge  and  687 
Ocean  Topography  Experiment  (TOPEX)/Poseidon  al¬ 
timetry  measurements. 


Table  la.  Statistics  of  water-level  elevation  (m)  and  phase 
(min)  for  PCTides  Yellow  Sea  predictions  compared  to 
observations. 


With  data  assimilation 

Water-level  elevation 

Phase 

RMSE 

<rMod 

oObs 

R 

SS 

Lag 

Location  1 

0.44 

0.52 

0.63 

0.72 

0.51 

0.95 

84 

Location  2 

0.43 

0.55 

0.58 

0.72 

0.47 

0.95 

84 

Location  3 

0.30 

0.60 

0.50 

0.87 

0.65 

0.94 

-48 

Location  4 

0.49 

1.06 

0.84 

0.89 

0.66 

0.94 

-36 

Average 

0.42 

0.68 

0.64 

0.80 

0.57 

0.95 

63* 

Without  data  assimilation 

RMSE 

aMod 

crObs 

R 

SS 

Rp 

Lag 

Location  1 

0.46 

0.63 

0.63 

0.73 

0.46 

0.95 

-84 

Location  2 

0.43 

0.63 

0.58 

0.75 

0.46 

0.95 

-84 

Location  3 

0.30 

0.50 

0.50 

0.82 

0.64 

0.94 

-60 

Location  4 

0.51 

0.87 

0.84 

0.82 

0.62 

0.94 

60 

Average 

0.43 

0.66 

0.64 

0.78 

0.55 

0.95 

72* 

The  asterisk  (*)  denotes  absolute  mean. 


PCTides  utilizes  tide  station  data  from  the  Interna¬ 
tional  Hydrographic  Office  (IHO)  database  (IHO 
1988)  for  model  validation  and/or  data  assimilation. 
The  IHO  database  in  PCTides  includes  tidal  data  for 
approximately  4500  tide  stations. 

PCTides  was  developed  to  generate  daily  to  weekly 
forecasts  of  water-level  elevations  and  tidally  driven 
currents  and  may  include  wind-driven  effects.  When 
wind  is  not  included,  the  forecasts  may  be  for  any  future 
time — next  week  or  next  year.  If  wind  is  included,  the 
forecast  corresponds  to  the  time  period  of  wind  forcing. 
For  most  navy  applications,  wind  and  surface  pressure 
are  obtained  from  the  U.S.  Navy's  Operational  Global 
Atmospheric  Prediction  System  (NOGAPS)  (Hogan 
and  Rosmond  1991)  or  the  Coupled  Ocean  Atmo¬ 
sphere  Mesoscale  Prediction  System  (COAMPS)  (Ho- 
dur  1997),  with  grid  resolutions  of  0.5°  and  0.2°,  respec¬ 
tively.  These  real-time  operational  wind  and  surface 
fields  are  available  via  the  navy's  Metcast  system 
(http://www.metnet.navy.mil/Metcast),  which  allows  a 
user  to  request  specific,  real-time,  atmospheric  forecast 
products.  A  second  source  of  wind  in  the  PCTides  sys¬ 
tem  is  provided  by  a  hurricane  model  (Holland  1980), 
which  produces  time-varying  hurricane  wind  fields 
based  on  a  prescribed  track,  minimum  surface  pressure, 
and  radius  of  maximum  winds.  PCTides  test  cases  using 
these  hurricane  winds  have  been  documented  for  Hur¬ 
ricane  Isabel  (Preller  et  al.  2005)  and  Hurricane  Ka¬ 
trina  (Blain  et  al.  2007).  A  third  source  of  wind  in  PC¬ 
Tides  is  the  manual  input  of  a  constant  value  for  wind 
speed  and  direction  over  the  entire  domain,  which  can 
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(a)  Station  1  (33.2700  N  125.2510  E) 


5  10  15  20  25  30 


September  1995 

Fig.  4.  (a)  Yellow  Sea  PCTides  depth-averaged  tidal  current  forecasts  (u  and  v  components;  ms1) 
with  data  assimilation  vs  the  depth-averaged  ADCP  observations  (black)  at  stations  1  and  2  for  the 
period  1-30  Sep  1995.  (b)  Yellow  Sea  PCTides  depth-averaged  tidal  current  forecasts  ( u  and  v  compo¬ 
nents;  m  s_1)  with  data  assimilation  vs  the  depth-averaged  ADCP  observations  (black)  at  stations  3  and 
4  for  the  period  1-30  Sep  1995. 


vary  temporally  but  not  spatially.  This  is  useful  in  the 
case  where  a  single,  local  wind  record  is  available.  It  is 
also  possible  to  run  PCTides  with  atmospheric  forcing 
from  other  meteorological  models;  however,  existing 
software  would  have  to  be  modified  to  process  the  new 
data  type. 


c.  PCTides  data  assimilation 

To  improve  the  simulation  of  tidally  forced  dynam¬ 
ics,  software  was  built  into  PCTides  to  (optionally) 
nudge  the  model  solution  to  tidal  heights  predicted  us¬ 
ing  data  from  IHO  tidal  stations  within  the  model  do- 
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(b)  Station  3  (35.9680  N  124.0330  E) 


5  10  15  20  25  30 


Fig.  4.  ( Continued ) 


main.  The  nudging  method  is  based  on  deriving  a  new 
solution  at  grid  points  near  each  tidal  station  from  a 
weighted  combination  of  the  model  solution  and  the 
station  observations.  The  weighting  function  is  calcu¬ 
lated  from  the  product  of  a  defined  weighting  constant 
and  the  following  equation: 

f3  =  cos(D  TT/dmax)  +  1,  (9) 

where  D  is  the  distance  from  the  grid  point  to  the  ob¬ 
servation  and  dmax  is  the  radius  of  influence,  which  is 


defined  as  40  km.  Figure  2  shows  an  example  of  I  HO 
locations  for  the  Yellow  Sea  region.  Global  1HO  tidal 
constituent  data  are  stored  in  ASCII  format,  and  new 
datasets  such  as  the  TOPEX/Poseidon  altimetry  data, 
formatted  appropriately,  can  easily  be  incorporated 
into  the  PCTides  tidal  dataset.  PCTides  has  been  tested 
with  and  without  data  assimilation.  In  some  locations 
with  open  coastlines,  the  model  performs  well  without 
data  assimilation.  However,  in  locations  with  more 
complex  coastal  geometry  and  bathymetry,  data  assimi- 
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Table  lb.  Statistics  of  tidal  u  and  u  components  (m  s  l)  and  phase  (min)  for  PCTides  Yellow  Sea  predictions  compared  to 

observations. 


With  data  assimilation 

U  (ms"1) 

U  phase 

V( 

ms1) 

V  phase 

RMSE 

o-Mod 

crObs 

R 

SS 

RP 

Lag  (min) 

RMSE 

o-Mod 

oObs 

R 

SS 

Rr 

Lag  (min) 

Location  1 

0.10 

0.19 

0.23 

0.90 

0.81 

0.97 

-48 

0.17 

0.30 

0.35 

0.88 

0.77 

0.98 

-60 

Location  2 

0.11 

0.18 

0.18 

0.81 

0.62 

0.95 

-72 

0.13 

0.29 

0.37 

0.95 

0.87 

0.98 

-36 

Location  3 

0.06 

0.13 

0.07 

0.99 

0.44 

0.99 

-12 

0.05 

0.12 

0.16 

0.98 

0.91 

0.99 

-12 

Location  4 

0.10 

0.28 

0.24 

0.94 

0.82 

0.99 

-36 

0.16 

0.22 

0.28 

0.81 

0.66 

0.98 

-72 

Average 

0.09 

0.20 

0.18 

0.91 

0.67 

0.98 

42* 

0.13 

0.23 

0.29 

0.91 

0.80 

0.98 

45* 

Without  data  assimilation 

U  ( ms  1 ) 

U  phase 

V(m  s  1 ) 

V  phase 

RMSE 

o-Mod 

crObs 

R 

SS 

RP 

Lag  (min) 

RMSE 

crMod 

oObs 

R 

SS 

Rr 

Lag  (min) 

Location  1 

0.16 

0.18 

0.23 

0.69 

0.46 

0.96 

-96 

0.19 

0.27 

0.35 

0.83 

0.69 

0.97 

-60 

Location  2 

0.20 

0.19 

0.18 

0.43 

-0.02 

0.96 

-132 

0.20 

0.28 

0.37 

0.85 

0.71 

0.97 

-60 

Location  3 

0.05 

0.06 

0.07 

0.76 

0.57 

0.97 

-84 

0.06 

0.17 

0.16 

0.95 

0.88 

0.98 

-36 

Location  4 

0.18 

0.18 

0.24 

0.66 

0.43 

0.96 

-96 

0.21 

0.25 

0.28 

0.70 

0.45 

0.96 

-96 

Average 

0.15 

0.15 

0.18 

0.64 

0.36 

0.96 

102* 

0.17 

0.24 

0.29 

0.83 

0.66 

0.97 

63* 

The  asterisk  (*)  denotes  absolute  mean. 


lation  often  improves  the  solution.  A  previous  study  by 
Blain  (1997)  showed  that  the  inclusion  of  coastal  tide 
data  can  improve  the  accuracy  of  both  coastal  and 
deep-water  tidal  prediction. 

d.  PCTides  execution  and  output 

A  typical  48-h  forecast  with  a  grid  resolution  of  5-10 
km  takes  approximately  2-10  min  to  run  using  an  Intel 
Core  2  CPU  6600  using  2.4  GHz  with  3.00  GB  of  RAM. 
PCTides  produces  2D  output  fields  of  water-level  el¬ 
evation  and  depth-averaged  ocean  currents  at  each 
model  grid  point.  In  addition,  water-level  elevation  and 
current  data  can  be  saved  at  user-selected  locations 
with  a  minimum  output  frequency  of  10  min.  The  out¬ 
put  can  be  viewed  as  an  ASCII  file  or  plotted  as  a  time 
series. 

3.  PCTides  validation 

PCTides  has  been  tested  in  over  15  regions  around 
the  globe  (Posey  et  al.  2006).  However,  only  a  small 
subset  of  these  comparisons  is  discussed  here.  Statisti¬ 
cal  evaluations  were  performed  on  the  PCTides  predic¬ 
tions  by  comparison  with  available  data,  including 
acoustic  Doppler  current  profilers  (ADCPs),  bottom- 
mounted  pressure  sensors,  and  National  Oceanic  and 
Atmospheric  Administration  (NOAA)  and  IHO  sta¬ 
tion  data.  The  error  statistics  used  include  root-mean- 
square  error  (RMSE),  standard  deviation  (cr),  correla¬ 
tion  coefficient  (/?),  and  skill  score  (SS)  (Kara  et  al. 
2003).  The  cross-correlation  function  ( Rp )  (Jenkins  and 
Watts  1969)  was  applied  to  the  tidal  phase  for  elevation 


and  velocity.  The  model’s  SS  is  computed  from  the  fol¬ 
lowing  equation: 

SS  =  R2  -  [R  -  (<Ty/(Tx)f  -  [Y  -  X)/(Txf .  (10) 

where  X  is  the  mean  of  the  data  values;  Y  is  the  mean 
of  the  model  values;  and  ox  and  crv  are  the  standard 
deviation  of  the  data  and  model  values,  respectively. 
An  SS  of  1.0  indicates  a  perfect  PCTides  prediction, 
and  a  negative  SS  indicates  that  PCTides  may  have 
normalized  amplitudes  larger  than  the  correlation  or 
large  biases  in  the  mean  (Murphy  and  Epstein  1989). 

To  illustrate  the  accuracy  of  the  system,  the  following 
section  describes  three  major  areas  where  PCTides 
have  been  validated  as  well  as  two  additional  cases  of 
PCTides  used  with  (Office  of  Naval  Research)  ONR- 
and  (North  Atlantic  Treaty  Organization)  NATO- 
funded  field  experiments. 

a.  Yellow  Sea  region 

Accurate  prediction  of  water-level  elevations  and 
currents  is  a  challenging  task  in  the  Yellow  Sea  because 
of  complex  bathymetry  and  coastal  geometry.  The  sea 
is  shallow  with  depths  ranging  from  90  m  along  the 
interior  trough  to  —20  m  within  50  km  of  the  coast. 
Tidal  amplitudes  vary  greatly  in  the  Yellow  Sea,  with 
high  ranges  (—9  m)  observed  near  Inchon.  During  the 
early  1990s,  NAVOCEANO  deployed  four  bottom- 
mounted  ADCP  moorings  in  the  Yellow  Sea  (Fig.  2).  A 
30-day  time  series  (1-30  September  1995)  of  velocity 
and  pressure  gauge  data  from  the  four  stations  was 
compared  against  two  PCTides  model  forecasts  without 
surface  forcing.  One  of  the  forecasts  assimilated  IHO 
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Fig.  5.  Bathymetry  (m)  of  the  Korean  Strait.  Squares  indicate  the  IHO  stations  and  circles  indicate 
the  location  of  bottom-mounted  ADCPs. 


data  and  the  other  did  not.  The  resolution  of  the  PC- 
Tides  grid  was  10  km. 

Figure  3  shows  a  comparison  of  water-level  elevation 
for  the  data  assimilative  case  versus  the  observations. 
Stations  1-3  show  tidal  ranges  of  ~2  m,  while  station  4 
(the  shallowest  location  and  closest  to  the  shore)  de¬ 
picts  a  range  of  almost  4  m.  Table  la  presents  the  water- 
level  elevation  statistics  calculated  using  the  model 
forecasts  with  and  without  data  assimilation  versus  the 
observations.  In  general,  the  statistics  for  the  water- 
level  elevation  using  data  assimilation  showed  a  slight 
improvement  ranging  from  2%  to  13%  over  the  results 
without  data  assimilation.  On  average,  the  RMSE  is 


lower,  while  the  R  and  SS  are  slightly  higher  in  the  data 
assimilation  case.  In  addition,  the  average  phase  lag, 
calculated  as  the  average  of  absolute  values,  is  9  min 
less  in  the  data  assimilation  case. 

The  depth-averaged  tidal  currents  from  the  PCTides 
forecasts  were  compared  against  the  depth-averaged 
ADCP  observations.  Figures  4a, b  show  comparisons  of 
the  depth-averaged  u-  and  u-velocity  components  of  the 
ADCPs  versus  PCTides  with  data  assimilation  for  the 
full  30-day  period.  The  phases  compare  well  quantita¬ 
tively  at  all  four  stations.  Stations  1  and  2  show  periods 
of  small  over-  or  underprediction  for  the  u  and  v  com¬ 
ponents  while  stations  3  and  4  show  overprediction  of 
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Station  8  (33.93  N,  128.34  E) 


Fig.  6.  Korean  Straits  PCTides  water-level  elevation  forecast  (red;  m)  with  data  assimilation  vs  the 
ADCP  (black)  at  stations  L2,  L4,  L6,  L8,  and  L10  for  the  period  1-11  Jun  1999. 


the  u  component  and  underprediction  of  the  v  compo¬ 
nent. 

Overall,  the  data  assimilative  case  showed  a  notice¬ 
able  improvement  in  R ,  SS,  and  lag  (see  Table  lb).  An 
examination  of  the  average  RMSE  of  the  individual  u 
and  v  components  showed  a  decrease  of  40%  and  24%, 
respectively,  with  the  inclusion  of  data.  The  effect  of 
data  assimilation  was  most  pronounced  for  the  average 


phase  lag  of  the  u  component,  with  a  reduction  of  1  h. 
The  v  component  showed  a  lag  time  reduced  from  63  to 
45  min.  Similar  to  the  water-level  elevations,  the  PC¬ 
Tides  tidal  current  forecasts  demonstrated  major  im¬ 
provements  when  data  assimilation  was  applied  to  the 
model.  Overall,  the  velocity  errors  and  phase  lags  are 
lower,  while  the  R  and  SS  are  higher  at  ADCP  locations 
3  and  4.  Water-level  elevation  statistics  for  the  model 
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show  the  same  trend  at  stations  3  and  4.  This  is  most 
likely  due  to  the  complex  bathymetry  data  near  the 
location  of  ADCPs  1  and  2,  which  is  not  resolved  by 
DBDB2  or  the  10-km  grid. 

In  general,  it  is  difficult  to  predict  tidal  conditions  in 
the  Yellow  Sea  because  of  its  complex  bathymetry, 
semienclosed  geometry,  and  large  tidal  signals.  Taking 
all  of  these  into  consideration,  the  globally  relocatable 
PCTides  model  showed  remarkable  skill  for  this  30-day 
period. 

b.  Korea/Tsushima  Strait  region 

The  Korea/Tsushima  Strait  lies  between  the  Repub¬ 
lic  of  Korea  and  Kyushu,  Japan,  and  is  a  shallow  area 
with  depths  of  —100  m.  Another  set  of  measurements 
used  in  the  evaluation  of  the  PCTides  forecasts  was 
located  in  this  region.  Figure  5  shows  the  positions  of  1 1 
bottom-mounted  ADCPs  deployed  for  11  months  dur¬ 
ing  1999-2000  as  part  of  NRL ’s  Dynamical  Linkage  of 
Asian  Marginal  Seas  (LINKS)  program  (Teague  et  al. 
20(H)).  The  moorings  were  divided  into  two  lines,  one 
on  the  Japan  Sea  side  of  the  strait  (L1-L5)  referred  to 
as  the  Northeast  ”  line  and  one  on  the  East  China  Sea/ 
Yellow  Sea  side  of  the  strait  (L6-L1 1)  referred  to  as  the 
Southwest  ”  line.  A  PCTides  domain  was  set  up  for  this 
region  with  4.7-km  resolution  and  run  for  the  period  1 
May-l  1  June  1999.  Two  test  cases  of  PCTides  forecasts 
were  run:  with  and  without  data  assimilation.  No  atmo¬ 
spheric  forcing  was  used  in  either  test  case. 

PCTides  forecast  water-level  elevations  and  currents 
were  compared  to  measurements  from  the  ADCPs. 
Teague  et  al.  (2006)  reported  larger  tidal  ranges  along 
the  southwest  line,  with  a  maximum  range  of  3.0  versus 
0.7  m  along  the  northeast  line.  Figure  6  shows  a  10-day 
comparison  of  observations  versus  PCTides-predicted 
water-level  elevations  with  data  assimilation  from  a 
sampling  of  the  moorings  (L2,  L4,  L6,  L8,  and  L10)  for 
the  period  M0  June  1999.  Smaller  water-level  eleva¬ 
tions  are  predicted  by  PCTides  at  L2  and  L4  (north¬ 
east)  and  larger  elevations  are  seen  at  L6,  L8,  and  L10 
(southwest)  in  agreement  with  the  observations.  The 
statistics  of  the  water-level  elevations  for  both  cases 
(with  and  without  data  assimilation)  are  shown  in  Table 
2a.  The  differences  between  the  two  regions  where  the 
ADCPs  were  deployed  are  evidenced  by  the  a  of  the 
data,  which  ranges  from  0.16  m  along  L1-L5  to  0.61  m 
along  L6-L11.  The  a  of  the  model  results  reproduces 
this  difference.  On  average,  the  tidal  phase  predictions 
with  data  assimilation  lagged  behind  the  observations 
by  22  min  along  the  northeast  line  and  by  18  min  along 
the  southwest  line,  compared  to  36  min  for  all  stations 
without  data  assimilation.  For  the  northeast  set  of 
moorings  (L1-L5)  the  data  assimilative  case  showed 


Table  2a.  Statistics  for  water-level  elevation  (m)  and  phase 
(min)  for  PCTides  for  Korean  Straits  predictions  compared  to 
observations. 

With  data  assimilation 

Water-level  elevation 

Phase 

RMSE  <rMod  oObs  R  SS 

RP  Lag 

Location  1 

0.04 

0.11 

0.13 

0.95 

0.89 

0.96 

-12 

Location  2 

0.05 

0.13 

0.15 

0.95 

0.90 

0.97 

-24 

Location  3 

0.05 

0.16 

0.16 

0.95 

0.91 

0.98 

-24 

Location  4 

0.05 

0.18 

0.18 

0.96 

0.91 

0.98 

-24 

Location  5 

0.07 

0.22 

0.20 

0.96 

0.88 

0.98 

-24 

Average 

0.05 

0.16 

0.16 

0.95 

0.90 

0.97 

22* 

Location  6 

0.14 

0.65 

0.67 

0.98 

0.96 

0.98 

-12 

Location  7 

0.13 

0.61 

0.63 

0.98 

0.96 

0.98 

-12 

Location  8 

0.13 

0.58 

0.60 

0.98 

0.95 

0.99 

-12 

Location  9 

0.13 

0.56 

0.58 

0.98 

0.95 

0.99 

-24 

Location  10 

0.13 

0.56 

0.58 

0.97 

0.95 

0.99 

-24 

Location  11 

0.14 

0.57 

0.58 

0.97 

0.94 

0.99 

-24 

Average 

0.13 

0.59 

0.61 

0.98 

0.95 

0.99 

18* 

Without  data  assimilation 

RMSE 

(rMod 

(jObs 

R 

SS 

RP 

Lag 

Location  1 

0.07 

0.17 

0.13 

0.93 

0.74 

0.98 

-36 

Location  2 

0.07 

0.18 

0.15 

0.93 

0.76 

0.98 

-36 

Location  3 

0.07 

0.20 

0.16 

0.94 

0.79 

0.98 

-36 

Location  4 

0.07 

0.21 

0.18 

0.95 

0.84 

0.98 

-36 

Location  5 

0.08 

0.23 

0.20 

0.95 

0.85 

0.98 

-36 

Average 

0.07 

0.20 

0.16 

0.94 

0.80 

0.98 

36* 

Location  6 

0.33 

0.86 

0.67 

0.94 

0.75 

0.99 

-36 

Location  7 

0.30 

0.80 

0.63 

0.94 

0.77 

0.99 

-36 

Location  8 

0.27 

0.75 

0.60 

0.95 

0.80 

0.99 

-36 

Location  9 

0.25 

0.72 

0.58 

0.95 

0.81 

0.99 

-36 

Location  10 

0.25 

0.71 

0.58 

0.95 

0.82 

0.99 

-36 

Location  11 

0.25 

0.72 

0.58 

0.95 

0.81 

0.99 

-36 

Average 

0.24 

0.76 

0.61 

0.95 

0.79 

0.99 

36* 

The  asterisk  (*)  denotes  absolute  mean. 


some  improvement  in  RMSE,  o\  R ,  SS,  and  phase  lag. 
However,  the  southwest  set  of  measurements  (L6-L1 1 ) 
showed  a  much  more  pronounced  improvement  with 
data  assimilation  indicated  by  a  46%  reduction  in 
RMSE,  22%  reduction  in  <r,  20%  improvement  in  SS, 
and  50%  reduction  in  phase  lag.  This  improvement  can 
be  attributed  in  part  to  the  larger  number  of  IHO  sta¬ 
tions  located  near  the  southwest  line  versus  the  north¬ 
east  line,  resulting  in  overall  better  agreement  with  the 
depth-averaged  ADCPs. 

Figures  7a, b  depict  comparisons  of  depth-averaged 
ADCP  currents  (both  u  and  u  components)  at  stations 
L3  and  L8  for  cases  with  and  without  data  assimilation. 
At  L3,  when  data  assimilation  is  included,  the  u  com¬ 
ponent  of  velocity  is  underpredicted  in  the  westward 
flow,  while  at  L8  the  data  assimilation  shows  an  overall 
better  comparison  with  data.  Table  2b  shows  the  statis¬ 
tical  comparison  of  the  observed  versus  the  predicted 
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(a)  Station  3  (Northeast)-  With  Data  Assimilation 


0.5 1 - 1 - r 

7 

t/> 


I 


>  -0.5' - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 — 

5  10  15  20  25  30  5  10 


Without  Data  Assimilation 


1 

3  _0.5 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 — 

5  10  15  20  25  30  5  10 


Fig.  7.  (a)  Korean  Straits  PCTides  depth-averaged  tidal  current  forecasts  ( u  and  v components;  ms-1) 
with  and  without  data  assimilation  vs  the  depth-averaged  ADCP  observations  (black)  at  station  3  for  the 
period  1  May-il  Jun  1999.  (b)  Korean  Straits  PCTides  depth-averaged  tidal  current  fore  casts  (u  and  v 
component;  ms-1)  with  and  without  data  assimilation  vs  the  depth-averaged  ADCP  observations 
(black)  at  station  8  for  the  period  1  May-11  Jun  1999. 


tidal  currents  generated  with  and  without  data  assimi¬ 
lation.  Book  et  al.  (2004)  found  inaccurate  bottom- 
depth  values  near  the  L6  mooring  when  using  the  Sung 
Kyun  Kwan  University  (Choi  1999)  dataset,  which  is 
included  in  the  PCTides  bathymetry.  As  a  result,  the  L6 
values  are  shown  but  not  used  in  the  averaged  statistics. 
The  average  statistics  for  the  u  and  v  components  show 


a  lower  RMSE  along  the  southwest  line  than  along  the 
northeast  line.  In  addition,  except  for  the  u  component 
in  the  data  assimilative  case,  the  phase  lag  is  smaller 
along  the  southwest  line.  Statistics  for  the  southwest 
line  are,  on  average,  improved  in  the  data  assimilation 
case.  Statistics  along  the  northeast  line  are  slightly  de¬ 
graded  in  the  data  assimilative  case. 
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(b)  Station  8  (Southwest)-  With  Data  Assimilation 


Without  Data  Assimilation 


As  an  oceanographic  constriction  or  chokepoint  be- 
tween  two  seas,  the  Korea  Strait  is  a  challenging  loca¬ 
tion  for  the  prediction  of  tides.  The  region,  one  of  com¬ 
plex  bathymetry  and  geometry,  has  been  observed  to 
have  two  distinct  tidal  conditions:  strong  tides  on  the 
East  China  Sea/Yellow  Sea  side  and  weak  tides  on  the 
Japan  Sea  side.  During  this  42-day  period,  the  PCTides 
forecast  water-level  elevation  was  able  to  duplicate 
these  conditions.  In  addition,  the  PCTides  currents, 
with  average  RMSE  ranging  from  0.03  to  0.07  ms-1 


and  average  phase  lags  ranging  from  10  to  67  min,  cap¬ 
tured  key  features  of  the  observed  currents  when  com¬ 
pared  to  the  two  lines  of  ADCP  observations. 

c.  Washington  State  region 

During  June  2005,  the  U.S.  Navy  conducted  an 
Autonomous  Undersea  Vehicle  Demonstration 
(AUVFest  2005),  held  at  the  U.S.  Navy  Undersea  War¬ 
fare  Center  in  Keyport,  Washington.  During  this  real¬ 
time  demonstration,  NAVOCEANO  performed  the  of- 
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Table  2b.  Statistics  of  tidal  u  and  v  components  (m  s 


)  and  phase  (min)  for  PCTides  Korean  Straits  predictions  compared  to 
observations. 


With  data  assimilation 


U  (m  s  l) 

U  phase 

K(ms  ') 

V  phase 

RMSE 

crMod 

crObs 

R 

SS 

Rn 

Lag  (min) 

RMSE 

crMod 

crObs 

R 

SS 

*p 

Lag  (min) 

Location  1 

0.05 

0.11 

0.12 

0.94 

0.86 

0.96 

-24 

0.07 

0.21 

0.16 

0.96 

0.78 

0.97 

-24 

Location  2 

0.06 

0.09 

0.12 

0.90 

0.76 

0.93 

-36 

0.06 

0.16 

0.15 

0.94 

0.86 

0.96 

-24 

Location  3 

0.06 

0.07 

0.12 

0.90 

0.71 

0.91 

-12 

0.06 

0.14 

0.13 

0.89 

0.75 

0.95 

-48 

Location  4 

0.07 

0.06 

0.11 

0.88 

0.64 

0.90 

36 

0.06 

0.13 

0.12 

0.92 

0.78 

0.97 

-48 

Location  5 

0.10 

0.09 

0.12 

0.69 

0.33 

0.92 

96 

0.06 

0.12 

0.13 

0.92 

0.78 

0.96 

-36 

Average 

0.07 

0.08 

0.12 

0.86 

0.66 

0.92 

41* 

0.06 

0.15 

0.14 

0.93 

0.79 

0.96 

36* 

Location  6** 

0.04 

0.18 

0.16 

0.98 

0.93 

0.98 

0 

0.04 

0.05 

0.05 

0.69 

0.33 

0.80 

96 

Location  7 

0.03 

0.14 

0.14 

0.98 

0.94 

0.99 

-12 

0.03 

0.04 

0.05 

0.83 

0.65 

0.84 

-36 

Location  8 

0.04 

0.13 

0.12 

0.98 

0.91 

0.99 

-12 

0.03 

0.06 

0.06 

0.89 

0.72 

0.95 

-60 

Location  9 

0.02 

0.11 

0.12 

0.98 

0.96 

0.99 

-24 

0.04 

0.07 

0.08 

0.86 

0.71 

0.98 

-84 

Location  10 

0.02 

0.10 

0.09 

0.99 

0.94 

0.99 

-12 

0.04 

0.08 

0.09 

0.90 

0.81 

0.99 

-72 

Location  11 

0.02 

0.12 

0.11 

0.98 

0.96 

0.98 

0 

0.06 

0.09 

0.09 

0.83 

0.63 

0.99 

-84 

Average 

0.03 

0.12 

0.12 

0.98 

0.94 

0.99 

10* 

0.04 

0.07 

0.07 

0.86 

0.70 

0.95 

67* 

Without  data  assimilation 

U  (m  s l) 

U  phase 

V  (m  s-1) 

V  phase 

RMSE 

crMod 

(xObs 

R 

SS 

K 

Lag  (min) 

RMSE 

o-Mod 

!  crObs 

R 

SS 

Lag  (min) 

Location  1 

0.06 

0.15 

0.12 

0.92 

0.74 

0.99 

-48 

0.10 

0.23 

0.16 

0.94 

0.63 

0.99 

-36 

Location  2 

0.06 

0.14 

0.12 

0.92 

0.75 

0.98 

-48 

0.06 

0.19 

0.15 

0.97 

0.86 

0.99 

-24 

Location  3 

0.05 

0.13 

0.12 

0.91 

0.79 

0.97 

-48 

0.06 

0.17 

0.13 

0.94 

0.76 

0.99 

-36 

Location  4 

0.05 

0.12 

0.11 

0.92 

0.84 

0.96 

-48 

0.05 

0.15 

0.12 

0.95 

0.81 

0.99 

-36 

Location  5 

0.06 

0.14 

0.12 

0.92 

0.75 

0.97 

-48 

0.05 

0.15 

0.13 

0.94 

0.84 

0.98 

-36 

Average 

0.06 

0.14 

0.12 

0.92 

0.77 

0.97 

48* 

0.06 

0.18 

0.14 

0.95 

0.78 

0.99 

34* 

Location  6** 

0.10 

0.24 

0.16 

0.96 

0.62 

0.97 

-12 

0.11 

0.10 

0.05 

-0.04 

-4.13 

0.40 

-192 

Location  7 

0.06 

0.19 

0.14 

0.97 

0.79 

0.98 

-24 

0.07 

0.08 

0.05 

0.53 

-0.96 

0.58 

72 

Location  8 

0.05 

0.16 

0.12 

0.97 

0.80 

0.98 

-24 

0.03 

0.08 

0.06 

0.92 

0.70 

0.93 

24 

Location  9 

0.03 

0.13 

0.12 

0.97 

0.93 

0.99 

-24 

0.02 

0.09 

0.08 

0.98 

0.92 

0.98 

-12 

Location  10 

0.02 

0.11 

0.09 

0.99 

0.94 

0.99 

-12 

0.02 

0.10 

0.09 

0.97 

0.93 

0.97 

0 

Location  11 

0.04 

0.12 

0.11 

0.95 

0.89 

0.99 

48 

0.03 

0.11 

0.09 

0.98 

0.91 

0.98 

-24 

Average 

0.04 

0.14 

0.12 

0.97 

0.87 

0.98 

24* 

0.03 

0.09 

0.07 

0.88 

0.50 

0.89 

26* 

The  asterisk  (*)  denotes  absolute  mean.  The  double  asterisk  (**)  denotes  that  the  station  is  not  used  in  average. 


ficial  operational  testing  (OPTEST)  of  the  PCTides  sys¬ 
tem.  At  the  end  of  the  evaluation  period,  the  model 
output  was  quantitatively  compared  to  the  ADCP  ob¬ 
servations  collected  during  the  exercise.  Prior  to  the 
OPTEST,  evaluation  criteria  defining  acceptable  fore¬ 
cast  errors  were  established  for  PCTides  to  meet  in 
order  to  be  declared  successful.  These  criteria  included 
1)  the  RMSE  of  PCTides-predicted  water-level  eleva¬ 
tion  versus  IHO  or  the  ADCP  amplitudes  relative  to 
mean  sea  level  must  be  less  than  0.30  m,  2)  the  RMSE 
of  PCTides-predicted  phase  versus  IHO  or  the  ADCP 
phase  must  be  less  than  30  min  (use  peak  amplitude 
times  only),  and  3)  the  RMSE  of  current  magnitude 
must  be  equal  to  or  less  than  0.5  ms-1. 

From  6  to  16  June  2005,  the  PCTides  system  fore¬ 
casted  water-level  elevations  and  depth-averaged  cur¬ 
rents  in  the  Washington  State  area.  The  first  grid  cov¬ 
ered  all  of  Puget  Sound  with  a  resolution  of  12  km  (Fig. 


8a).  A  higher-resolution  nested  grid  was  run  daily  with 
a  resolution  of  1.5  km  (Fig.  8b).  Three  test  cases  were 
run:  1)  with  data  assimilation,  2)  without  data  assimila¬ 
tion,  and  3)  with  assimilation  and  atmospheric  wind/ 
pressure.  FES99  tidal  constituents  were  applied  as  lat¬ 
eral  boundary  conditions  to  the  PCTides  large  outer 
grid  (Fig.  8a)  while  the  nest  received  boundary  con¬ 
ditions  from  the  outer  grid  (12  km).  The  bathymetry 
utilized  for  AUVFest  2005  included  NRLs  DBDB2 
and  a  high-resolution  3-s  bathymetry  provided  by 
NAVOCEANO.  PCTides  software  was  used  to  blend 
these  datasets  to  produce  an  improved  bathymetry  and 
provide  more  accurate  model  forecasts. 

Similar  to  the  Yellow  Sea  and  Korea  Straits  test 
cases,  PCTides  simulations  were  performed  with  and 
without  data  assimilation.  Additional  tests  were  run  to 
investigate  the  effect  of  wind  forcing  with  differing 
spatial  resolution.  Winds  used  during  these  tests  were 
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Fig.  8.  (a)  Host  grid  resolution  of  12  km.  (b)  Nestl  grid  resolution  of  1.5  km  with  IHO  stations 
shown  and  (c)  zoomed  area  of  the  Hood  Canal  with  ADCP  shown. 


from  the  COAMPS  model  available  at  resolutions  of 
27  km  (a  standard  operational  resolution)  and  a  higher 
resolution  of  5  km  (typically  available  only  by  special 
request).  Our  testing  showed  that  the  most  accu¬ 
rate  forecasts  were  made  utilizing  winds  from  the 
higher-resolution  (5  km)  test  case.  Subsequently, 
only  the  higher-resolution  wind  case  results  will  be 
shown. 

The  validation  data  used  in  this  study  are  from  an 
ADCP  deployed  in  the  Hood  Canal  by  NAVOCEANO 
during  the  period  7-18  June  2005.  The  ADCP  data  in¬ 


dicated  that  the  tidal  range  in  this  region  was  greater 
than  4  m.  Figure  9  depicts  a  water-level  comparison  of 
observations  versus  the  model  forecasts  for  three  test 
cases:  1 )  without  data  assimilation  or  wind,  2)  with  data 
assimilation  and  no  wind,  and  3)  with  data  assimilation 
and  wind.  Table  3a  presents  water-level  elevation  sta¬ 
tistics  for  these  three  test  cases.  The  wind  test  case 
showed  the  best  overall  skill  with  an  RMSE  of  0.16  m, 
which  is  impressive  considering  the  tidal  range  in  this 
region  was  greater  than  4  m.  The  Rp  for  the  winds  case 
was  0.99  with  a  lag  of  20  min.  While  the  data  assimila- 


770 


JOURNAL  OF  ATMOSPHERIC  AND  OCEANIC  TECHNOLOGY 


Volume  25 


Hood  Canal  -  Without  Winds  and  Without  Data  Assimilation 


June  2005 


Fig.  9.  Water-level  elevation  (m)  comparisons  of  Hood  Canal  ADCP  (black)  vs  PCTides  forecasts 
(red)  using  (a)  no  wind  and  no  data  assimilation,  (b)  no  wind  and  data  assimilation,  and  (c)  wind  and  data 
assimilation.  The  plot  is  valid  6-17  Jun  2005.  PCTides  grid  resolution  is  1.5  km.  The  wind  and  data 
assimilation  case  started  1  day  later  than  the  other  two  cases. 


tion  cases  compared  well  against  the  observations,  the 
study  without  data  assimilation  did  poorly.  The  RMSE 
was  ~1  m  while  the  R  and  SS  were  low.  In  these  test 
cases,  the  PCTides  output  frequency  was  chosen  to  be 
20  min;  as  a  result,  a  smaller  phase  lag  could  not  be 
resolved. 

PCTides  also  forecast  tidal  current  magnitudes  dur¬ 


ing  AUVFest  2005  as  part  of  the  operational  evaluation 
criteria  discussed  previously.  The  tidal  current  magni¬ 
tude  comparison  between  the  ADCP  and  the  PCTides 
forecast  from  the  1.5-km  grid  using  1 )  no  wind  and  data 
assimilation  and  2)  COAMPS  wind  with  data  assimila¬ 
tion  is  shown  in  Fig.  10.  The  statistics  for  both  cases  are 
very  similar.  The  RMSE  with  the  inclusion  of  wind  was 
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Table  3a.  Statistics  of  water-level  elevation  (m)  and  phase 
(min)  for  PCTides  Hood  Canal  predictions  compared  to 
observations. 


With  data  assimilation 

Water-level  elevation 

Phase 

RMSE 

trMod  (rObs  R  SS 

Lag 

Hood  Canal  with  wind 

0.16 

0.95  1.05  0.99  0.97 

0.99 

20 

Hood  Canal  without  wind 

0.20 

1.00  1.05  0.98  0.96 

0.99 

20 

Without  data  assimilation 

RMSE 

crMod  (jObs  R  SS 

Rp 

Lag 

Hood  Canal  without  wind 

0.95 

1.28  1.05  0.69  0.20 

0.87 

96 

0.05  m  s"1  and  without  was  0.06  ms-1.  The  R  for  both 
cases  was  0.84. 

During  AUVFest  2005,  PCTides  forecasts  were  used 
to  provide  an  optimum  dive  window  in  two  recovery 
missions  (a  buoy  and  a  bottom  crawler).  Unfortunately, 
because  of  poor  visibility,  the  buoy  was  not  rescued.  A 
bottom  crawler,  which  had  become  stuck  in  the  mud, 
was  successfully  recovered.  Afterward,  the  feedback 
from  the  divers  was  that  the  PCTides  model  forecasts 
were  dead  on.  "  In  addition  to  this  feedback,  the  sta¬ 
tistical  comparison  calculated  from  the  PCTides  fore¬ 
casts  versus  the  observations  was  used  to  determine 
that  PCTides  met  all  the  OPTEST  criteria  and  was  de¬ 
clared  operational  by  NAVOCEANO. 

Although  not  part  of  the  PCTides  AUVFest  evalua¬ 
tion  criteria,  we  examined  the  u  and  v  components  of 
the  depth-averaged  current  from  the  model  versus  the 
depth-averaged  ADCP  currents  (similar  to  the  com¬ 
parisons  discussed  earlier  for  the  Yellow  Sea  and  Ko¬ 
rean  straits).  The  Hood  Canal  presents  a  challenging 
environment  for  tidal  prediction  with  a  narrow  channel 
and  limited  IHO  stations  for  data  assimilation  (see  Fig. 
8b).  Figure  11  depicts  comparisons  of  the  depth- 
averaged  ADCP  currents  (both  u  and  v  components)  at 
the  Hood  Canal  location  for  all  three  PCTides  cases. 
Table  3b  shows  the  statistical  comparison  of  the  u  and 
u  components  of  the  PCTides  depth-averaged  current 
versus  the  observation  of  the  Hood  Canal  ADCP.  Both 
the  u  and  v  components  from  the  data  assimilation  and 
wind  case  had  high  R  and  SS  values  with  a  phase  lag  of 
0  min.  The  case  with  data  assimilation  but  without  wind 
had  higher  RSME  and  lower  SS  than  the  previous  case. 
Similar  to  the  statistics  for  the  water-level  elevation,  the 
velocity  statistics  for  the  case  without  data  assimilation 
or  wind  showed  the  poorest  results.  These  results  show 
the  importance  of  data  assimilation  in  constricted  areas 
as  well  as  wind  forcing. 


d.  Other  exercises  (ASIAEX  and  MREA04) 

During  the  spring  of  2001,  the  Asian  Seas  Interna¬ 
tional  Acoustics  Experiment  (ASIAEX)  took  place  in 
the  South  and  East  China  Seas.  To  study  soliton  gen¬ 
eration  near  the  Luzon  Strait,  Ramp  et  al.  (2004)  ran 
PCTides  and  compared  water-level  elevations  to  obser¬ 
vations  at  two  locations.  The  model  was  run  for  an  area 
from  15°-26°N,  115°-125°E  and  assimilated  data  from 
54  IHO  stations.  In  this  earlier  version  of  the  model, 
FES95. 1/2.1  (Shum  et  al.  1997)  was  used  to  provide 
lateral  boundary  conditions.  The  authors  found  excel¬ 
lent  agreement  in  water-level  elevations  during  the  pe¬ 
riod  22  April-19  May  2001  with  an  RMSE  of  0.15  m. 
The  authors  also  compared  the  PCTides  results  against 
the  Oregon  State  University  TOPEX/Poseidon  global 
inverse  solution  (Egbert  et  al.  1994)  and  found  that 
PCTides  produced  ‘Virtually  identical  results  ”  in  the 
examination  of  water-level  elevations  for  the  area  of 
interest. 

PCTides  was  the  tidal  prediction  component  of 
NRL's  wave,  tide,  and  current  prediction  system  used 
during  the  NATO  Maritime  Rapid  Environmental  As¬ 
sessment  (MREA04)  trial  held  in  Portuguese  coastal 
waters  near  Pinheiro  da  Cruz  during  April  2004  (Fig. 
12).  The  PCTides  grid  was  set  up  in  this  area  with  a 
resolution  of  0.05°  and  provided  daily  48-h  forecasts  of 
water-level  elevations.  Wind  and  sea  level  pressure 
from  COAMPS  (27  km)  were  used  as  atmospheric  forc¬ 
ing  for  the  daily  forecasts.  A  total  of  28  IHO  stations 
were  assimilated  into  the  model  solutions. 

Waves  and  nearshore  currents  were  measured  during 
MREA04  by  deploying  a  Nortek  Vector  acoustic  velo- 
cimeter  (N6)  in  the  surf  zone.  The  sensor  package  in¬ 
cluded  a  pressure  gauge  placed  at  low  tide  in  1.5  m  of 
water.  Figure  13  shows  a  scatterplot  of  water-level  el¬ 
evations  from  PCTides  versus  water  observations  at 
Pinheiro  da  Cruz,  an  area  where  surf  predictions  for 
amphibious  landings  were  being  generated  during 
MREA04.  The  mean  water  depth  from  the  N6  data  was 
subtracted  from  the  measured  water  depth  to  deter¬ 
mine  the  water-level  elevation.  Since  the  N6  data  were 
recorded  at  irregular  intervals  (60  or  120  min),  the  PC¬ 
Tides  water-level  elevations  were  temporally  interpo¬ 
lated  to  match  the  N6  observation  times.  The  statistics 
showed  excellent  agreement  with  an  RMSE  in  ampli¬ 
tude  of  0.07  m,  phase  errors  less  than  15  min,  and  an  R 
of  approximately  1.0. 

4.  Summary  and  conclusions 

A  global  tide/surge  forecast  system,  PCTides,  has 
been  developed  by  NRL  for  the  U.S.  Navy  to  predict 
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Hood  Canal  ADCP  (47.785-122.732) 
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a) 


June  2005 


June  2005 


Fig.  10.  Tidal  current  magnitude  (m)  comparisons  of  Hood  Canal  ADCP  (black)  vs  PCTides  forecasts 
(red)  using  (a)  no  wind  and  data  assimilation  and  (b)  COAMPS  wind  and  data  assimilation.  The  plot  is 
valid  6-47  Jun  2005.  PCTides  grid  resolution  is  1.5  km.  The  wind  and  data  assimilati  on  case  started  1  day 
later  than  the  other  case. 


water-level  elevation  and  depth-averaged  ocean  cur¬ 
rents.  One  of  the  challenges  for  the  navy  is  to  have  the 
ability  to  predict  water-level  elevation  and  currents 
for  essentially  any  region  worldwide.  Tidal  predic¬ 
tion  in  shallow,  semienclosed  seas  and  constricted  re¬ 
gions  such  as  straits  and  channels  is  often  more  diffi¬ 
cult  than  prediction  in  open  coastal  regions.  To  prove 


the  global  robustness  of  the  PCTides  system,  it  has 
been  evaluated  in  many  different  locations  by  a  number 
of  users.  This  study  has  focused  on  three  of  the  more 
challenging  locations  for  tidal  prediction:  the  shallow, 
semienclosed  Yellow  Sea,  the  Korea  Strait,  and  the 
Hood  Canal.  PCTides  exhibited  skill  in  prediction 
of  both  water-level  elevation  and  tidal  currents  at 
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Hood  Canal  -  Without  Winds  and  Without  Data  Assimilation 


>6  7  8  9  10  11  12  13  14  15  16  17  18 


Winds  and  Data  Assimilation 
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June  2005 


Fig.  1 1.  Hood  Canal  PCTides  depth-averaged  tidal  current  forecasts  (w  and  v components;  m  s'1)  with 
(a)  no  wind  and  no  data  assimilation,  (b)  no  wind  and  data  assimilation,  and  (c)  wind  and  data  assim¬ 
ilation  vs  the  depth-averaged  ADCP  observation  (black)  at  Hood  Canal.  The  plot  is  valid  6-17  Jun  2(X)5. 
PCTides  grid  resolution  is  1.5  km.  The  wind  and  data  assimilation  case  started  1  day  later  than  the  other 

two  cases. 


all  three  locations.  The  Hood  Canal  case  was  used  as 
the  U.S.  Navy’s  operational  testing  of  PCTides.  The 
system  easily  met  the  evaluation  criteria  defined  by  the 
U.S.  Navy  in  this  location.  It  was  also  shown  that  in 
these  three  locations,  the  PCTides  tidal  forecasts  were 
improved  by  assimilating  IHO  tidal  station  data.  This 
holds  true  for  most  locations  in  which  PCTides  has 


been  tested.  This  was  particularly  evidenced  in  the  wa¬ 
ter-level  elevation  statistics  for  the  southwest  transect 
in  the  Korean  Straits  and  Hood  Canal  where  the  as¬ 
similation  of  IHO  station  data  showed  a  marked  im¬ 
provement  in  model  skill.  Although  this  paper  focused 
primarily  on  the  evaluation  of  PCTides  both  with  and 
without  data  assimilation,  the  Hood  Canal  test  case 
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Table  3b.  Statistics  of  tidal  u  and  u  components  (m  s  *)  and  phase  (min)  for  PCTides  Hood  Canal  predictions  compared  to 


observations. 

With  data  assimilation 

U  ( ms  ') 

U  phase 

K(ms  ') 

V  phase 

RMSE 

crMod 

crObs 

R 

SS 

K 

Lag  (min)  RMSE  crMod  oObs 

R 

SS 

RP 

Lag  (mm) 

Hood  Canal  with  wind 

0.04 

0.09 

0.07 

0.92 

0.71 

0.93 

0  0.04 

0.11  0.11 

0.94 

0.89 

0.94 

0 

Hood  Canal  without  wind 

0.07 

0.14 

0.12 

0.93 

0.18 

0.93 

20  0.05 
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showed  that  the  inclusion  of  wind  forcing  added  to  the 
forecast  skill. 

This  study  has  illustrated  the  importance  of  accurate 
bathymetry  as  well  as  the  value  of  data  assimilation  in 
tidal  prediction.  Although  the  globally  relocatable  PC- 
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Fig.  12.  Bathymetry  (m)  of  MREA04  exercise  off  the  coast  of 
Portugal.  The  star  is  the  Nortek  acoustic  velocimeter  (N6)  near 
Pinheiro  da  Cruz. 


Tides  system  has  been  proven  to  be  a  skillful  tool  for 
tidal  prediction,  continued  feedback  from  users  of  the 
system  and  subsequent  model  enhancements  will  fur¬ 
ther  improve  the  system’s  accuracy. 
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